APS/123-QED 



New measurement of neutron capture resonances of 209 Bi 

C. Domingo-Pardo 1 *, U. Abbondanno 2 , G. Aerts 3 , H. Alvarez-Pol 4 , F. Alvarez- Velarde 5 , S. Andriamonje 3 , 
J. Andrzejewski 6 , P. Assimakopoulos 7 , L. Audouin 8 , G. Badurek 9 , P. Baumann 10 , F. Becvaf 11 , E. Berthoumieux 3 , 
F. Calvino 12 , D. Cano-Ott 5 , R. Capote 13 ' 14 , A. Carrillo de Albornoz 15 , P. Cennini 16 , V. Chepel 17 , 
E. Chiaveri 16 , N. Colonna 18 , G. Cortes 12 , A. Couture 19 , J. Cox 19 , M. Dahlfors 16 , S. David 20 , I. Dillman 8 , 
R. Dolfini 21 , W. Dridi 3 , I. Duran 4 , C. Eleftheriadis 22 , M. Embid-Segura 5 , L. Ferrant 20 , A. Ferrari 16 , 
R. Ferreira-Marques 17 , L. Fitzpatrick 16 , H. Frais-Koelbl 13 , K. Fujii 2 , W. Furman 23 , R. Gallino 24 , I. Goncalves 17 , 
E. Gonzalez- Romero 5 , A. Goverdovski 25 , F. Gramegna 26 , E. Gricsmayer 13 , C. Guerrero 5 , F. Gunsing 3 , B. Haas 27 , 
R. Haight 28 , M. Heil 8 , A. Herrera-Martinez 16 , M. Igashira 29 , S. Isaev 20 , E. Jericha 9 , Y. Kadi 16 , F. Kappeler 8 , 

D. Karamanis 7 , D. Karadimos 7 , M. Kerveno, 10 , V. Ketlerov 25 ' 16 , P. Koehler 30 , V. Konovalov 23 ' 16 , 
E. Kossionides 31 , M. Krticka 11 , C. Lamboudis 7 , H. Leeb 9 , A. Lindote 17 , I. Lopes 17 , M. Lozano 14 , S. Lukic 10 , 

J. Marganiec 6 , L. Marques 15 , S. Marrone 18 , P. Mastinu 26 , A. Mengoni 13 ' 16 , P.M. Milazzo 2 , C. Moreau 2 , 
M. Mosconi 8 , F. Neves 17 , H. Oberhummer 9 , M. Oshima 32 , S. O'Brien 19 , J. Pancin 3 , C. Papachristodoulou 7 , 
C. Papadopoulos 33 , C. Paradela 4 , N. Patronis 7 , A. Pavlik 34 , P. Pavlopoulos 35 , L. Perrot 3 , R. Plag 8 , A. Plompen 36 , 
A. Plukis 3 , A. Poch 12 , C. Pretel 12 , J. Quesada 14 , T. Rauscher 37 , R. Reifarth 28 , M. Rosetti 38 , C. Rubbia 21 , 
G. Rudolf 10 , P. Rullhusen 36 , J. Salgado 15 , L. Sarchiapone 16 , I. Savvidis 22 , C. Stephan 20 , G. Tagliente 18 , 
J.L. Tain 1 , L. Tassan-Got 20 , L. Tavora 15 , R. Terlizzi 18 , G. Vannini 39 , P. Vaz 15 , A. Ventura 38 , D. Villamarin 5 , 
M. C. Vincente 5 , V. Vlachoudis 16 , R. Vlastou 33 , F. Voss 8 , S. Walter 8 , H. Wendler 16 , M. Wiescher 19 , K. Wisshak 8 



The n_TOF Collaboration 

1 Instituto de Fisica Corpuscular, CSIC-Universidad de Valencia, Spain 
2 Istituto Nazionale di Fisica Nucleare, Trieste, Italy 
3 CEA/Saclay - DSM, Gif-sur-Yvette, France 
4 Universidade de Santiago de Compostela, Spain 
5 Centro de Investigaciones Energeticas Medioambientales y Technologicas, Madrid, Spain 

6 University of Lodz, Lodz, Poland 
7 University of Ioannina, Greece 
8 Forschungszentrum Karlsruhe GmbH (FZK), Institut fur Kernphysik, Germany 
9 Atominstitut der Osterreichischen Universitdten, Technische Universitdt Wien, Austria 
10 Centre National de la Recherche Scientifique/IN2P3 - IReS, Strasbourg, France 
11 Charles University, Prague, Czech Republic 
12 Universitat Politecnica de Catalunya, Barcelona, Spain 
13 International Atomic Energy Agency, NAPC-Nuclear Data Section, Vienna, Austria 
14 Universidad de Sevilla, Spain 
15 Instituto Tecnologico e Nuclear (ITN), Lisbon, Portugal 
16 CERN, Geneva, Switzerland 
17 LIP - Coimbra & Departamento de Fisica da Universidade de Coimbra, Portugal 
18 Istituto Nazionale di Fisica Nucleare, Bari, Italy 
19 University of Notre Dame, Notre Dame, USA 
20 Centre National de la Recherche Scientifique/IN2P3 - IPN, Orsay, France 
21 Universita degli Studi Pavia, Pavia, Italy 
22 Aristotle University of Thessaloniki, Greece 
23 Joint Institute for Nuclear Research, Frank Laboratory of Neutron Physics, Dubna, Russia 
24 Dipartimento di Fisica, Universita di Torino and Sezione INFN di Torino, Italy 
25 Institute of Physics and Power Engineering, Kaluga region, Obninsk, Russia 
26 Istituto Nazionale di Fisica Nucleare(INFN), Laboratori Nazionali di Legnaro, Italy 
27 Centre National de la Recherche Scientifique/IN2P3 - CENBG, Bordeaux, France 
28 Los Alamos National Laboratory, New Mexico, USA 
29 Tokyo Institute of Technology, Tokyo, Japan 
30 Oak Ridge National Laboratory, Physics Division, Oak Ridge, USA 
31 NCSR, Athens, Greece 
32 Japan Atomic Energy Research Institute, Tokai-mura, Japan 
33 National Technical University of Athens, Greece 
34 Institut fur Isotopenforschung und Kernphysik, Universitat Wien, Austria 
35 Pole Universitaire Leonard de Vinci, Paris La Defense, France 
36 CEC-JRC-IRMM, Geel, Belgium ' 
37 Department of Physics and Astronomy - University of Basel, Basel, Switzerland 

3H ENEA, Bologna, Italy 
39 Dipartimento di Fisica, Universita di Bologna, and Sezione INFN di Bologna, Italy^ 

(Dated: February 6, 2008) 



2 



The neutron capture cross section of Bi has been measured at the CERN n_TOF facility 
by employing the pulse height weighting technique. Improvements over previous measurements 
are mainly due to an optimized detection system, which led to a practically negligible neutron 
sensitivity. Additional experimental sources of systematic error, like the electronic threshold in 
the detectors, summing of 7-rays, internal electron conversion, and the isomeric state in bismuth, 
have been taken into account. 7- Ray absorption effects inside the sample have been corrected 
by employing a non-polynomial weighting function. Since 209 Bi is the last stable isotope in the 
reaction path of the stellar s-process, the Maxwellian averaged capture cross section is important 
for the recycling of the reaction flow by a decays. In the relevant stellar range of thermal energies 
between kT = 5 keV and 8 keV our new capture rate is about 16% higher than the presently 
accepted value used for nucleosynthesis calculations. At this low temperature an important part of 
the heavy Pb-Bi isotopes are supposed to be synthesized by the s-process in the He shells of low 
mass, thermally pulsing asymptotic giant branch stars. With the improved set of cross sections we 
obtain an s-process fraction of 19±3% of the solar bismuth abundance, resulting in an r-process 
residual of 81±3%. The present (n, 7) cross section measurement is also of relevance for the design 
of accelerator driven systems based on a liquid metal Pb/Bi spallation target. 

PACS numbers: 25.40.Lw,27.80.+w,97.10.Cv 



I. INTRODUCTION 

209 Bi is the end point isotope of the s-process path. 
Its observed abundance has been finally understood in 
terms of the main s-process component, operating in the 
He-shell of low mass, low metallicity, thermally pulsing 
asymptotic giant branch (agb) stars pi Q • 

However, the complex production pattern at the end 
of the s-process path still hides important information, 
that can be unraveled by means of accurate neutron cap- 
ture measurements. Capture on bismuth leads to the 
ground state of 210 Bi, which /3-decays to 210 Po that is a- 
unstable, producing 206 Pb (see Fig. QJ. However, given 
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FIG. 1: (Color online) The s-process reaction network ter- 
minates at 209 Bi. Shaded boxes designate stable isotopes. 
Decays from the a-unstable 210 > 211 Po isotopes are represented 
with dashed lines. Quantities below the isotopic symbols show 
natural abundances or half lifes. 



the 210 Po half-life of 138 days it can also act as a branch- 
ing point by capturing another neutron and enhancing 
the abundance of 207 Pb instead. Another contribution 
to 207 Pb may be due to a long-lived isomer in 210 Bi at 
271 keV. Neutron capture on this isomer leads to 211 Bi 
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with an cv-decay half-life of 2.17 min. The two branch- 
ings, at 210 Po and 210 Bi, depend strongly on the stellar 
conditions of neutron density and temperature. Given an 
appropriate stellar model and accurate neutron capture 
cross sections, the s abundances of the Pb/Bi isotopes 
can be more reliably determined, thus providing a bet- 
ter decomposition of the respective r-process abundances 
and of the radiogenic contributions due to the Th/U a 
decays. The latter information represents a viable con- 
straint on the Th/U abundances calculation and its use 
as cosmochronometers 0, • 

For the discussion of these astrophysical aspects, the 
present status of the 209 Bi capture data is rather un- 
satisfactory Previous experiments exhibit significant 
discrepancies suggesting the existence of systematic un- 
certainties, which must be clearly reduced for a more 
quantitative assessment of the Pb/Bi abundances. Two 
capture resonance studies on 209 Bi have been reported so 
far, one at ORNL @ and the other at GELINA 0,@. In 
the former measurement, two strong s-wave resonances 
at low neutron energy could not be measured because 
of a low neutron energy cutoff at 2.6 keV. In addition, 
substantial corrections were needed for some of the res- 
onances in order to account for the neutron sensitivity 
of the experimental setup employed. The GELINA mea- 
surement was also affected by significant corrections for 
neutron sensitivity. 

The 209 Bi (n, 7) cross section has also a very practical 
application. Thanks to its properties of chemical inert- 
ness, high boiling point, low neutron moderation, and 
large scattering cross section, an eutectic mixture of lead 
and bismuth is presently considered as a very appropri- 
ate material for the spallation target and as the coolant 
for accelerator driven systems (ADS). 

In a first sensitivity study of the neutron cross sec- 
tions of the main materials used in this type of hybrid 
reactor Q , severe discrepancies where found between dif- 
ferent evaluated nuclear data files. As a consequence, 
a list of isotopes was tagged for high priority measure- 
ments. Initially, most of the requested cross sections cor- 
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responded to the fuel and cladding materials, and were fo- 
cused on the elastic, inelastic and (n, xn) channels. How- 
ever, a recent extension of the same study includes 
additional requests for the capture reactions in the spal- 
lation/coolant materials lead and bismuth. Comparing 
two evaluated data files, ENDF/B-VI.8 and JENDL 3.3, 
discrepancies in the neutron balance of 12.5% were iden- 
tified to result from the uncertainties of the (n, 7) cross 
section of 209 Bi alone. Hence, the precise knowledge of 
the neutron capture cross sections of the lead and bis- 
muth isotopes turned out to be of key relevance for the 
design of an ADS suited for the transmutation of radioac- 
tive residues and for energy production. 

Apart from its contribution to the neutron balance, 
neutron capture on 209 Bi affects the radiotoxicity as well. 
Build-up of the a-emitters 210 Bi m (t 1/2 = 3xl0 6 yr) and 
210 Po (fi/2 = 138 d) contributes to the long and short 
term radiotoxicity of the spallation target and of the 
coolant. These hazards can only be evaluated on the 
basis of a precise in, 7) cross section of bismuth and with 
accurate data concerning the partial cross sections lead- 
ing to 210 Bi m and 210 Po. 



II. CROSS SECTION MEASUREMENT 

The neutron energy differential capture cross section 
of 209 Bi has been measured at the n_TOF installation of 
CERN using the time of flight (TOF) technique. Since 
this facility has been described in detail elsewhere [Tlj . 
only the main features will be summarized in section Al Al 
The major motivation for this measurement was to over- 
come the systematic uncertainty related to the neutron 
scattering background, which arises in measurements on 
isotopes with a dominant scattering channel. Sample 
scattered neutrons are easily captured in the materials 
surrounding the detector or in the sensitive detector vol- 
ume itself, producing 7-rays, which appear as a prompt, 
TOF-dependent background. This effect becomes partic- 
ularly crucial for the broad s-wave resonances in bismuth 
at low neutron energies between 0.8 keV and 12 keV, re- 
sulting in corrections as large as 50% in previous mea- 
surements |6j. 

For this reason the reaction yield of this isotope has 
been measured with a set of total energy detectors con- 
sisting of optimized CqQq liquid scintillation detectors 
( section [II B|) . This setup turned out to be of pivotal im- 
portance for the efficient suppression of the background 
from scattered neutrons. Other sample-related correc- 
tions were significantly reduced by choosing a compara- 
bly small bismuth sample 20 mm in diameter and 6.1 mm 
in thickness. The cross section is determined relative to 
the 4.9 eV 198 Au resonance by employing the saturated 
resonance technique . For this purpose a gold sample 
1 mm in thickness and 20 mm in diameter was measured 
in regular intervals. A thorough analysis of the capture 
data, including extensive Monte Carlo simulations of the 
necessary corrections to the experimental data, allowed 



us to determine the final cross section with a systematic 
uncertainty of better than 3%. 

A. The n_TOF installation 

At n_TOF, neutrons are produced via spallation reac- 
tions induced by 20 GeV protons in a lead target. Pulses 
of (3-7) x 10 12 protons with a duration of 6 ns (rms) hit 
the Pb-block with a typical repetition rate of 0.4 Hz. 
This very low duty cycle combined with the very high in- 
tensity in the beam pulses makes n_TOF an ideal facility 
for TOF measurements on radioactive samples. The low 
repetition rate has allowed to implement an advanced 
data acquisition system with zero dead time based on 
Flash-amplitude-to-digital converters |l3j . 

The number of neutrons arriving at the sample are 
monitored by means of a 200 ng/cm 2 thick 6 Li-foil, 
centered in the beam. The particles emitted in the 
6 Li(n, a) 3 H reaction are registered with four silicon de- 
tectors surrounding the 6 Li-foil outside of the beam. The 
neutron monitor [14J is mounted in an evacuated carbon 
fiber chamber 3 m upstream of the samples. The sam- 
ples, which are positioned at a flight path of 185 m, are 
fixed on a sample changer, which is also made of carbon 
fiber. 

The neutron intensity has been determined by means 
of two independent measurements, performed with the 
Si-monitor described above and with a calibrated fission 
chamber |l5j |. Two measurements were carried out with 
the latter detector, employing foils of 235 U and 238 U. Cal- 
ibration of the absolute yield in (n, 7) measurements via 
the saturated resonance method requires actually only 
the shape of the neutron flux, which can be directly de- 
termined from the Si-monitor and fission chamber mea- 
surements with an uncertainty below ±2% |l6l . 

The beam profile at the sample position has been de- 
termined to exhibit an approximate Gaussian shape with 
a = 7 mm, slightly off centered (A^ = 1.5 mm ) [17|- The 
samples used in this measurement cover about 60% of the 
neutron beam. 

The excellent neutron energy resolution of the n_TOF 
facility is due to the long flight path of 185 m and to the 
short proton pulse width. The resolution function was 
determined by means of Monte Carlo simulations an d 
has been validated experimentally by measuring narrow 
p- and ei-wave capture resonances on a sample of nat Fe. 



B. Experimental setup 

A general view of the experimental setup is shown in 
Fig. The detection system used with the pulse height 
weighting technique (PHWT) 0] had been optimized 
with respect to the sensitivity for scattered neutrons, the 
most crucial source of background in (n, 7) studies on 
neutron magic nuclei, which are characterized by partic- 
ularly large scattering to capture ratios. Scattered neu- 
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trons may be subsequently captured in the detectors or in 
nearby materials, where they produce 7-rays, which are 
likely to be registered with high probability, thus mim- 
icking true capture events in the sample. 

In an effort to reduce the neutron sensitivity as much 
as possible, special CqDq liquid scintillation detectors pcj 
have been developed consisting of materials with low neu- 
tron capture cross sections. The main features of these 
detectors are thin-walled scintillator cells made of carbon 
fiber, which are directly glued onto the photomultiplier 
tubes, thus eliminating dead materials, e.g. the common 
quartz window of the scintillator cells. To avoid the addi- 
tional material of a support, the detectors were hanging 
on thin cords fixed at the ceiling. Similarly, the sample 
changer was made of carbon fiber as was the frame of the 
sample-holder. The samples were mounted on thin Kap- 
ton foils, which were glued on the frame that was much 
larger than the diameter of the neutron beam. In this 
way, the neutron sensitivity in the critical energy range 
from 1 keV to 100 keV could be reduced by factors of 
three to ten pcj . 

The two detectors were placed at an angle of ^125° 
with respect to the sample in order to minimize the an- 
gular distribution effects of the primary capture 7 rays. 
Also the background due to in-beam 7-rays scattered in 
the sample was considerably reduced in this configuration 



Eq.^are the better justified the better conditions (i) and 
(ii) are fulfilled. 

The validity of Eq. ^ can be challenged by several ex- 
perimental problems though. In particular, condition (i) 
is violated if more than one 7-ray of the cascade is regis- 
tered in the same detector. Moreover, the product as well 
as the sum over the m 7-rays of the cascade is always in- 
complete because of the unavoidable loss of 7-rays due to 
electronic threshold applied to the detector signals, due 
to converted transitions, and due to partial population 
of the isomeric state in 210 Bi at 271 keV. 

For all these effects, which may influence the validity 
of Eq. ^ appropriate corrections have to be determined 
as is shown in the following subsections. 

As far as condition (ii) is concerned, the proportional 
increase of the efficiency with 7-ray energy is enforced 
by an appropriate modification of the detector energy 
response distribution R(E). This is achieved by appli- 
cation of a pulse height dependent weighting factor Wi , 
such that the weighted sum of the response for a 7-ray 
j , Rij , becomes proportional to its energy E 1 j , 

J2WiRij =aE yj . (2) 

i 

The second approximation in Eq. ^ depends directly on 
the accuracy of the calculated weights Wi , which can be 
tested as discussed below. 



III. DATA ANALYSIS 

In this section, the PHWT principle will be reviewed 
together with the related systematic uncertainties of rele- 
vance for this type of capture measurements, followed by 
the determination of the weighting factors. We describe 
the calculation of accurate yield corrections and the anal- 
ysis procedure to determine resonance parameters. 

A. Pulse height weighting technique 
and systematic uncertainties 

The PHWT is based on two conditions, (i) that the 
detector efficiency £ 7 < 1 so that at most only one 7- 
ray per capture cascade is registered, and (ii) that e 7 is 
proportional to the energy of the registered 7-ray, e 7 « 
otE-f. Under these two conditions, the efficiency e c for 
detecting a capture event, i.e. a cascade composed of m 
7-rays, becomes proportional to the sum energy E c of 
that cascade. In this case one obtains 



£ c = 1 - Y[(l -e 7 j) « X) s, w ~ aEc > W 

which is a constant independent of the actual de- 
excitation pattern of the nucleus produced in the capture 
reaction. It is worth noting that the approximations in 



B. Weighting factors 

It has been shown [22I |23| that realistic response func- 
tions Rij for mono-energetic 7-rays with energies E 1 j can 
be determined by means of Monte Carlo simulations for 
any particular setup for (n, 7) measurements. A set of 
expressions of the form of Eq. [21 can then be used to de- 
rive the weighting factors Wi. 




FIG. 2: (Color online) The geometry of the experimental 
setup used in the Geant4 simulation of the weighting func- 
tions. 

The detailed geometry and chemical composition of 
the experimental setup has been implemented in the 
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Monte Carlo simulation using Geant4 [24| (see Fig. 01 . 
The composition of the carbon fiber, which is by 
far the dominant structural material, was determined 
by an RBS analysis [H, yielding C/O/N/Ca/Br = 
2.0/0.2/0.16/0.012/0.016. Also, the capture 7-rays were 
carefully traced, assuming a radial distribution for the 
emission probability inside the sample according to the 
neutron beam profile described in section lil Al The depth 
distribution was included as well since it changed signif- 
icantly from a practically uniform to a surface peaked 
shape between weak and strong resonances, respectively. 

The weighting function (WF) of the gold calibration 
sample was obtained using the conventional approxima- 
tion by a polynomial function of order k = 4, Wi = 

J2k=t ak ^i ■ The varues °f the coefficients Ofc were de- 
rived from a least squares minimization, 



min 




Because of the higher 7-ray absorption in the 6.1 mm 
thick bismuth sample the proportionality condition (ii) 
could not be satisfied with the polynomial WF approach. 
In this case an accurate, pointwise WF was obtained by 
means of a linear regularization method for solving the 
inverse problem of Eq.|2 16] . In this way the uncertainty 
connected with the polynomial WF has been reduced to 
0.3%, an improvement by an order of magnitude. 

The systematic uncertainty introduced by the WF was 
determined as described in Ref. 23], basically by per- 
forming Monte Carlo simulations of the capture cascades 
at a certain resonance. The de-excitation pattern is mod- 
eled by combining the experimentally known levels at low 
excitation energy with a statistical model of the nucleus 
to complete the cascade up to the capture energy. These 
simulated capture events can be used in order to estimate 
the uncertainty of the calculated WF. 

In order to illustrate the performance of the cascade 
event generator for the samples measured in this exper- 
iment, the simulated pulse height spectra Rf for the 
802 eV resonance in 210 Bi and for the 4.9 eV resonance 
in 198 Au are compared in Fig. [3] with the corresponding 
experimental spectra measured at n_TOF. 




E deD (MeV) E dep (MeV) 



FIG. 3: Comparison of experimental (dashed) and simulated 
(solid) pulse height spectra of 210 Bi (left) and 198 Au (right) 
illustrating the quality of the simulations. 

The spectra based on the simulated cascades show in 
general a good agreement with the experimental spec- 



tra. The small differences at higher energy can be as- 
cribed to the employed statistical model and to a lesser 
extent to the uncertainties in the calibration of the GqDq 
detectors. In principle, these differences could be mini- 
mized by modifying the relative strength of Ml and El 
transitions in the statistical part of the cascades and/or 
the instrumental resolution. However, the effect on the 
weighted sum ^ ■ WiRf and, therefore, on the estimated 
uncertainty becomes absolutely negligible in both cases. 

In summary, it was found that the simulated pulse 
height spectra can be employed for the determination of 
the uncertainty of the WF. By computing the weighted 
sum of each capture cascade with the calculated weight- 
ing functions, we obtain values, which deviate by less 
than 0.3% from the corresponding capture energies. This 
means that condition (ii) is fulfilled within 99.7%, leading 
to a practically negligible uncertainty due to the calcu- 
lated WF. 



C. Yield correction factors 

The corrections for the electronic threshold, 7- 
summing, internal conversion, and the effect of the iso- 
meric state (See Sec. IIII A|) have been estimated by de- 
tailed Monte Carlo simulations of the capture events in 
a similar manner as in Ref. |23|. While internal conver- 
sion and summing of cascade 7-rays were found to have 
a minor effect of less than 2% on the capture yield of 
the bismuth sample, more significant corrections are re- 
quired for the electronic threshold and for the isomeric 
state in 210 Bi. The isomeric state is included in the sim- 
ulation by assigning a null branching ratio to the level of 
271 keV, so that the de-excitation cascade ends at that 
level whenever it is reached. 

The yield correction factor including all these effects 
can be calculated as 



where Rf is the total response of the detection sys- 
tem for n simulated capture events. For a threshold t of 
200 keV, the yield correction factors calculated with this 
procedure for the bismuth sample are listed in the third 
column of Table [I] 

Since the simulated spectra for resonances with same 
spin and parity are very similar, i.e. almost independent 
of the resonance energy, the correction factors could be 
classified according to spin and parity of the resonance 
(first column in Table QJ. 

The population of the isomeric state calculated with 
our statistical model of the nucleus is shown in the sec- 
ond column of TableQ] For comparison, and to illustrate 
the possible uncertainty of the calculation, the correction 
factor of the hypothetic case, where the isomer is pop- 
ulated with 100% probability is shown in the last col- 



6 



TABLE I: Yield correction factors for the bismuth sample 
calculated for a threshold of 200 keV. 



J w isom. P Opul.(%) ft.s.ce.m j-t,s,ce,V 



3+ 6 1.124(2) 1.103(2) 

4" 8 1.131(2) 1.118(2) 

4+ 9 1.143(2) 1.116(2) 

5~ 23 1.137(2) 1.109(2) 

5+ 22 1.142(2) 1.119(2) 

6+ 44 1.134(2) 1.114(2) 



umn. The differences were found to range between 1% 
and 2.5%. 

As mentioned before, the cross section is determined 
relative to the 4.9 eV resonance in 197 Au. Therefore, 
the capture yield measured with the gold sample had to 
be corrected for experimental effects of threshold, 7-ray 
summing, and internal conversion. The corresponding 
correction factor for a threshold of 200 keV is f^f e = 
1.046(2). 

Hence, the final correction for each resonance is the 
ratio between the corresponding factor yM. ce < m from Ta- 
ble U and fj^f e - This ratio oscillates between 5% and 
7%. It can be concluded, that the threshold, 7-summing 
and internal conversion effects do not cancel out by mea- 
suring with respect to a reference sample. Neglecting 
these corrections could lead to systematic deviations in 
the yield determination of >5%. 



where e 7 is the detection probability for a 7-ray of a 
given energy. For an energy E 1 = 600 keV, £„/e 7 has 
been calculated by means of Monte Carlo simulations for 
comparison with the experimentally determined value re- 
ported in Ref. 20]. The second factor, e 7 /e c , has been 
determined by a detailed Monte Carlo simulation of the 
present experimental setup. The efficiency e 7 was ob- 
tained for 7-rays of 600 keV, whereas e c was determined 
using simulated capture cascades (see sections fill Bl and 
1111 C|) . Finally, r„/r 7 is to be determined by an iterative 
procedure of correcting and fitting the value of T 7 of the 
corresponding resonance. 

In case of the s-wave resonance at 12.1 keV, which has 
the largest r„/r 7 of all analyzed bismuth resonances, the 
corresponding ratios were found to be e 7 /e c = 0.446 (this 
work) and (e n /e 7 )i2.i keV = 1-5841 x 10~ 5 [2(j. Because 
of the very low neutron sensitivity the required correc- 
tion for the 12.1 keV resonance of 3.6±0.7% is relatively 
small, and practically negligible for the rest of the s-wave 
resonances in bismuth. 

The 20% uncertainty of this correction results from 
neutron captures in the quartz window of the photomul- 
tiplier tube due to the uncertain silicon cross sections, 
which contribute a 10% uncertainty to e„/e 7 , and from 
the resonance parameters r„ and r 7 , which exhibit un- 
certainties of 8% and 10%, respectively. 



E. R-matrix analysis 

The experimental yield 



D. Neutron sensitivity 



For a certain resonance at energy E with capture and 
neutron widths, T 7 and T n , the probability that a signal 
in the C§T>e scintillator is caused by a neutron scattered 
in the sample and eventually captured in the detector, is 
given by 



pns — I zlL \ I Z2L 



(5) 



where e n denotes the probability to detect a 7-ray pro- 
duced by sample-scattered neutrons in or near the detec- 
tor and e c the probability to register a 7-ray from a true 
capture cascade. This requires a correction of the exper- 
imental resonance yield by 



1 + P" 



Eq. [SJcan be rewritten as, 



pns _ 



(7) 



-y-exp jt.s,ce.in jns jsat 



N" 



N n E, 



(8) 



was determined by the weighted net count rate {N w ), 
the effective binding energy E Cl the integrated neutron 
flux N n (obtained from the Si-monitor and from the 
shape of the neutron flux, sec. Ill A| . and the corrections 
discussed before. The factor f sat corresponds to the ab- 
solute normalization via the analysis of the saturated res- 
onance at 4.9 eV in 198 Au. 

The yield has been analyzed with the multilevel R- 
matrix code SAMMY (2(|. Using an iterative procedure 
based on Bayes' theorem, Y exp is fitted with a function 
of the type, 



Y ex ? = Yf(E ,r n ,r^) + B, 



(9) 



where Y* corresponds to a parameterization as a func- 
(6) tion of the resonance parameters E , T n and T 7 accord- 
ing to the Reich-Moore formalism. For all partial waves a 
channel radius of 9.6792 fm was used. While adopting the 
well known T n values from transmission measurements, 
the resonance energies E and capture widths T 7 have 



been fitted to the experimental data. The background 
below the resonance is described by the constant term 
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B. In the fits with SAMMY also sample effects (sin- 
gle and double neutron scattering inside the sample) as 
well as thermal broadening are taken into account. The 
resolution function of the facility (sec. Ill A]l has been con- 
sidered in the fits using the RPI parameterization. 

The present results are illustrated in Fig. 0] at the 
example of the first two s-wave resonances in 210 Bi. 
The comparison with the capture yields calculated with 
the corresponding resonance parameters quoted in the 
ENDF/B-VI.8 evaluation demonstrates the improvement 
due to the reduced systematic uncertainties of the n_TOF 
data. 



0.01 




E„ (eV) E n (keV) 



FIG. 4: (Color online) R-matrix analysis of the first two 
resonances in bismuth. The dashed line corresponds to the 
yield calculated with the resonance parameters quoted by the 
ENDF/B-VI.8 evaluation. 



TABLE II: Resonance parameters" and radiative kernels 6 for 
209 Bi. 



E (eV) 


I 


J 


T n (meV) 


T 7 (meV) 


gr 7 r„/r (meV) 


801.6(1) 





5 


4309(145) 


33.3(12) 


18.2(6) 


2323.8(6) 





4 


17888(333) 


26.8(17) 


12.0(8) 


3350.83(4) 


1 


5 


87(9) 


18.2(3) 


9.5(2) 


4458.74(2) 


1 


5 


173(13) 


23.2(22) 


11.3(11) 


5114.0(3) 





5 


5640(270) 


65(2) 


35.3(11) 


6288.59(2) 


1 


4 


116(18) 


17.0(17) 


6.7(7) 


6525.0(3) 


1 


3 


957(100) 


25.3(14) 


8.6(5) 


9016.8(4) 


1 


6 


408(77) 


21.1(14) 


13.0(9) 


9159.20(7) 


1 


5 


259(45) 


21.4(21) 


10.9(11) 


9718.910(1) 


1 


4 


104(22) 


74(7) 


19.5(21) 


9767.2(3) 


1 


3 


900(114) 


90(8) 


28.7(26) 


12098 










65 (4) c 


15649.8(1.0) 


1 


5 


1000 


47(4) 


20.2(17) 


17440.0(1.3) 


1 


6 


1538(300) 


32(3) 


20.4(18) 


17839.5(9) 


1 


5 


464(181) 


43(4) 


21.7(20) 


20870 


1 


5 


954(227) 


34.4(33) 


18.3(17) 


21050 


1 


4 


7444(778) 


33(3) 


14.8(13) 


22286.0(9) 


1 


5 


181(91) 


33.6(32) 


15.1(15) 


23149.1(1.3) 


1 


6 


208(154) 


25.3(25) 


14.7(15) 



"Angular orbital momenta, I, resonance spins J. and 
neutron widths, V n , are mainly from Refs. [27l |2^|. 
b Uncertainties are given as 18.2(6) = 18.2±0.6. 
c This area corresponds to the sum of the areas of the broad 
s-wave resonance at the indicated energy, plus two p-wave 
resonances at 12.092 and 12.285 keV. 



IV. 



RESULTS 



In total, 21 resonances were identified in the 209 Bi(n, 7) 
data in the energy range from 800 eV up to 23150 eV. The 
respective capture widths and areas are listed in Table ITU 

Beyond 23 keV further resonances could not be ob- 
served due to in-beam 7-ray background |2ll | , which was 
the major limitation in this experiment. However, this 
background did not affect the initial aim of improving the 
systematic uncertainty of the broad s-wave resonances 
below 100 keV. 

The capture areas determined at n_TOF are compared 
in Fig. [3 with the previous measurements performed at 
ORNL 6] and GELINA 0. 

For most of the resonances shown in Fig. there is 
a fair agreement with the ORNL data. However, sig- 
nificantly higher values are reported in the latter mea- 
surement for the two resonances at 3.3 keV and 5.1 keV, 
and a lower cross section for the resonance at 9.77 keV. 
The GELINA data show agreement with n_TOF in sev- 
eral cases, but there is an average trend of lower cross 
sections. The largest differences in absolute values are 
for the resonances at 5.1 keV, 9.77 keV and the group 
at 12.1 keV. For these three resonances, the ORNL data 
shows opposite trend for the 5.1 keV resonance, agrees 
with GELINA for the 9.77 keV and agrees with the 
n_TOF data for the 12.1 keV group. 

The general systematics of better average agreement 
with the ORNL data and lower average values of 
the GELINA data, might be related to the weighting 
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FIG. 5: (Color online) Comparison of the radiative capture 
kernels derived from the measurements at n_TOF, GELINA 
and ORNL. The capture area shown at 2.3 keV for ORNL 
corresponds to an estimation reported in Ref . . The capture 
area shown at ~12keV, corresponds to the sum of the three 
resonances in that region. 



functions used in the respective experiments, although 
threshold correction effects (see Sec. IIII C|) could play a 
role as well. Systematic deviations due to the neutron 
sensitivity corrections applied to some particular reso- 
nances can not be excluded also. 
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A. Discussion of uncertainties 

The measuring technique employed in this work and 
the data analysis procedures described in Sec. II I II and 
IIII CI have been validated experimentally at n_TOF [22| • 
According to this study, the various contributions to the 
systematic uncertainty of the present data are about 2%. 
In addition, the energy dependence of the neutron flux 
has been found to exhibit an uncertainty of 2% as well, 
resulting in a total systematic uncertainty of better than 
3%. 

In this discussion, the neutron sensitivity correction 
has to be considered separately since it has been applied 
only to one of the measured resonances. Although the 
uncertainty of the correction was 20%, the contribution 
to the systematic uncertainty of the capture kernel was 
only 0.7% in this case. 

V. IMPLICATIONS FOR THE s- ABUNDANCES 
IN THE PB-BI REGION 

The slightly higher cross sections for s-wave resonances 
measured at n_TOF with respect to the values obtained 
at GELINA may affect the synthesis of the Pb/Bi. This 
holds preferentially at lower stellar temperature, since 
the neutron sensitivity correction becomes less relevant 
above 100 keV, where the relative contribution of the 
neutron scattering channel decreases. 

Qualitatively, the consequence of a higher 209 Bi cross 
section would reduce the survival probability and, hence, 
the s-process abundance of this isotope, but this could be 
compensated to some extent through an enhanced pro- 
duction through a-recycling via 206 > 207 pb. 

In the calculation of the Maxwellian averaged cross sec- 
tions (MACS), the n_TOF resonances have been comple- 
mented between 25 keV and 31 keV with resonances from 
Ref. , and in the interval from 31 keV up to 60 keV with 
capture areas from Ref. [27j ■ The additional resonances 
represent 2% and 7% of the MACS at kT = 5 keV and 
kT = 8 keV, respectively. At lower stellar temperature, 
comparison of the results in Table lllll with the MACS 
obtained at GELINA @ shows that the new values are 
16% larger. 



TABLE III: Maxwellian averaged cross sections of Bi com- 
pared to recent previous data. 



kT 


This work 


Mutti et al.[8j 


(keV) 


(mb) 


(mb) 


5 


13.05(84) 


11.25(58) 


8 


8.62(54) 


7.48(44) 


20 


3.41(50) 


3.34(42) 


25 


2.89(50) 


2.85(45) 



Beyond kT = 20 keV, the MACS is strongly influenced 
by the average capture cross section above 80 keV, which 
has been determined experimentally in Ref. For this 



reason, the differences obtained in the strong s-wave reso- 
nances at lower neutron energy do not affect substantially 
the MACS at 20-25 keV. 

The MACS given in Table lTTTl refer to the lower temper- 
ature regime of the common stellar s-process site associ- 
ated with thermally pulsing low mass AGB stars [2^ . Ac- 
cording to the Galactic chemical evolution (GCE) model 
described in Refs. [J, the s-process abundances of 
208 Pb and 209 Bi are essentially produced in stars of low 
metallicity. In this model about 95% of the neutron expo- 
sure is due to the 13 C(a,n) 16 reaction, which operates 
during the interpulse phase between He shell flashes at 
temperatures around ~ 10 8 K, corresponding to a ther- 
mal energy of kT w 8 keV. 

The additional neutron irradiation provided by the 
22 Ne(a, n) 25 Mg reaction at the higher thermal energy of 
kT = 23 keV during the He shell flash is rather weak. 
Hence, the abundances of isotopes with small cross sec- 
tions are still dominated by what has been produced at 
the lower temperature of the interpulse phase. Therefore, 
the low temperature part of the MACS is important for 
this stellar model. With respect to the solar abundance 
tables of Anders and Grevesse [3(|, the GCE model Q 
provides s-process abundances of 206 Pb, 207 Pb and 209 Bi 
at the epoch of the solar system formation of 62%, 79%, 
and 19%, respectively. 

According to the sensitivity study reported in Ref. |f| , 
the present result (see table IIIl|) is not expected to affect 
the abundances of 206 ' 207 Pb by a-recycling from A = 210. 
This was confirmed by a model calculation for a ther- 
mally pulsing AGB star with M = 3Mq and a metallicity 
[Fe/H] = —1.3, which yields a negligible difference for 
the abundances of 206 Pb and 207 Pb. The s-process abun- 
dance of 209 Bi itself was found to decrease only slightly 
from 19% to 18.7%. 

To evaluate the uncertainty on the solar 209 Bi s-process 
fraction, we have to consider first the uncertainty in 
the cross section of this isotope in the two temperature 
regimes of the stellar model. The main contribution due 
to the uncertainty of the bismuth cross section is domi- 
nated by the conditions during the He shell flash. At the 
higher temperature of the 22 Ne(a,n) 25 Mg source (kT = 
23 keV) the reaction flow via 208 Pb(rt, 7) 209 Bi is strongly 
enhanced over the situation during the 13 C(a,n) 16 
phase. This reflects the increase of the 208 Pb cross section 
with stellar temperature j^l). The net effect of the un- 
certainty on the 209 Bi cross section to its solar s-process 
abundance is found to be 6%. 

In the evaluation of the final s-process abundance of 
bismuth, three aditional uncertainties have to be consid- 
ered: (i) the effect of the uncertainty in the 208 Pb cross 
section, which directly affects the s production of 209 Bi 
by 6.5% (see Table VII in Ref. Q), (ii) the 7-8% uncer- 
tainty of the solar bismuth abundance (Refs. [3(], 15^ 1. 
and (iii) a further 10% for the uncertainties related to 
the s-process model and for modeling GCE 0. In sum- 
mary the s-process contribution to solar bismuth is ob- 
tained as 19±3%, corresponding to an r-process residual 
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of 81±3%. 

This result is in agreement with the r-process calcula- 
tion of Ref. [13 using the waiting point approximation 
and improved mass formulae, which yield an r-process 
contribution between 71% and 90%. A similar result is 
reported in Ref. where the calculations give 77% and 
92% r-process contribution depending on the initial seed 
composition. 

VI. ESTIMATION OF THE THERMAL 
CAPTURE CROSS SECTION 

The thermal neutron capture cross section of 209 Bi can 
be expressed by the sum of the tails of all Breit-Wigner 
resonances, 

4.099 x 10 6 (^j £<p^- (10) 

The result obtained with the resonance parameters 
from the present measurement (see Table ITU, is in agree- 
ment with that obtained from the GELINA measure- 
ment, as shown in Table IIVI However, both values 
are around 40% smaller than the accepted value of 
33.8(5) mb [27|], which has been measured with the pile 
oscillator method. The direct capture process cannot ac- 
count for this discrepancy. In fact, a preliminary estimate 
of this component leads to a negligible contribution to the 
thermal cross section value j34(. 

TABLE IV: Comparison of the thermal neutron capture cross 
section from different sources. 

Ref. al] 1 (mb) 

In-pile measurement [2JJ 33.8(5) 

From resonance parameters 

This work 23.6(9) 
GELINA [7] 24.6(9) 

ENDF 32.51 
JENDL 32.51 



The main improvement with respect to previous mea- 
surements is due to an optimized detection setup, by 
which the crucial neutron sensitivity could be consider- 
ably reduced. All remaining sources of systematic uncer- 
tainties have been thoroughly treated by detailed Monte 
Carlo simulations. 

Resonance energies (E ), widths (r 7 ), and capture ar- 
eas have been determined by an R-matrix analysis cov- 
ering the energy range from 0.8 keV to 23 keV. 

The results show larger capture areas for the s-wave 
resonances, yielding a 16% enhancement of the stellar 
neutron capture rate at thermal energies between 5 keV 
and 8 keV compared with recent data 0- 

The new cross section results for 209 Bi, combined with 
recent improvements in the cross section of 208 Pb p|, 
yield a solar s-process abundance of 19(3)% for bismuth. 
The resulting r-process residual of 81(3)% represents a 
reliable constraint for r-process calculations ; 33j. 
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The resonance parameters of the ENDF and JENDL 
libraries have been adjusted to reproduce the in-pile a^ h 
measurement by means of an expression like Eg. 1101 With 
the improved, present data for the s-wave resonances 
this inconsistency could be removed by introducing sub- 
threshold resonances in the evaluated data files. 



VII. SUMMARY 



At the CERN n_TOF facility the time of flight method 
has been employed with the pulse height weighting tech- 
nique in order to determine the neutron capture cross 
section of bismuth in the resolved resonance region. 
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